Thermalization process in bare and dressed coordinate approaches 
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We consider a particle in the approximation ol a harmonic oscillator, coupled linearly to a field 
modeling an environment. The field is described by an infinite set of harmonic oscillators, and the 
system (particle-field) is considered in a cavity at thermal equilibrium. We employ the notions of 
bare and dressed coordinates to study the time evolution of the occupation number. With dressed 
coordinates no renormalization procedure is required, leading directly to a finite result. In particular, 
for a large time, the occupation number of the particle becomes independent of its initial value. So 
we have a Markovian process, describing the particle thermalization with the environment. 

A thermalization process occurs in some cases for a system of material particles coupled to an environment, in the 
sense that after an infinitely long time, the matter particles loose the memory of their initial states. This study is, in 
general, not easy from a theoretical point of view, due to the complex non-linear character of the interactions between 
the matter particles and the environment. To get over these difficulties, linearized models have been adopted. An 
account on the subject of the evolution of quantum systems on general grounds can be found in [l|, @, U 0, H, @] . 
Besides, the main analytical method used to treat these systems at zero or finite temperature is, except for a few special 
cases, perturbation theory. In this framework, the perturbative approach is carried out by means of the introduction 
of bare, non-interacting objects (fields, to which are associated bare quanta), the interaction being introduced order 
by order in powers of the coupling constant. 

In spite of the remarkable achievements of the perturbative methods, however, there are situations where they 
cannot be employed, or are of little use. These cases have led to attempts to improve non-perturbative analytical 
methods, in particular, where strong effective couplings are involved. Among these trials there are methods that 
perform resummations of perturbative series, even if they are divergent, which amounts in some cases to extending 
the weak-couplin g reg ime to a strong-coupling domain. One of these methods is the Borel resummation of perturbative 
series @,i i El ETUI- 

In this paper we follow a different non-perturbative approach: we investigate a simplified linear version of a particle- 
field or particle-environment system, where the particle, taken in the harmonic approximation, is coupled to the 
reservoir, modeled by independent harmonic oscillators @, S We will employ, in particular, dressed states and 
renormalized coordinates. Using this method non-perturbative treatments can be considered for both weak and strong 
couplings. A linear model permits a better understanding of the need for non-perturbative analytical treatments of 
coupled systems, which is the basic problem underlying the idea of a dressed quantum mechanical system. Of course, 
the use of such an approach to a realistic non-linear system is an extremely hard task, while the linear model provides 
a good compromise between physical reality and mathematical reliability. The whole system is supposed to reside 
inside a spherical cavity of radius R in thermal equilibrium at temperature T — /3 . In other words, we consider 
the spatially regularized theory (finite R) at finite temperature. The free space case is obtained by suppressing the 
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regulator, (R_ 
space, see [13| 



oo). For a detailed comparison between this procedure and the one considering an a priori unbounded 



I. THE MODEL 

Let us start by considering a particle approximated by a harmonic oscillator, having bare frequency ujq, linearly 
coupled to a set of N other harmonic oscillators, with frequencies u>k, k = 1,2, . . . , N. The Hamiltonian for such a 
system is written in the form, 
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leading to the following equations of motion, 
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q\ +ujfqi = Ciq (t). 



(2) 
(3) 



In the limit N — > oo, we recover our case of the particle coupled to the environment, after redefining divergent 
quantities, in a manner analogous to mass renormalization in field theories. A Hamiltonian of the type ([1]) has been 
largely used in the literature, in particular to study the quantum Brownian motion with the path-integral formalism [l|, 
HJ . It has also been employed to investigate the linear coupling of a particle to the scalar potential [H, Q [l5|, [JJ|, \Vq . 
The Hamiltonian |T]) is transformed to principal axis by means of a point transformation, 
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JV 



r=0 r=0 

fi = (0,{fc}), fe=l,2,...,iV; r = 0,...N, 



(4) 



performed by an orthonormal matrix T — (t*). The subscripts ji = and fi = k refer respectively to the particle and 
the harmonic modes of the reservoir and r refers to the normal modes. In terms of normal momenta and coordinates, 
the transformed Hamiltonian reads 



1 N 



(5) 



r=0 



where the f2 r 's are the normal frequencies corresponding to the collective stable oscillation modes of the coupled system. 
Using the coordinate transformation ((J) in the equations of motion and explicitly making use of the normalization of 
the matrix (t r ), £^ =0 (^) 2 = 1, we get 
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with the condition 



(6) 



(7) 



We take c k = T](u>k) u , where rj is a constant independent of fc. In this case the environment is classified according to 
u > 1, u = 1, or u < 1, respectively as supraohmic, ohmic or subohmic. This terminology has been used in studies of 
the quantum Brownian motion and of dissipative systems [1, [H; H, [|| • For a subohmic environment the sum in Eq. Q 
is convergent in the limit N — ► oo and the frequency lj q is well defined. For ohmic and supraohmic environments, this 
sum diverges for N — > oo. This makes the equation meaningless, unless a renormalization procedure is implemented. 
From now on we restrict ourselves to an ohmic system. In this case, Eq. ([7]) is written in the form 
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(8) 
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where we have defined the counterterm 

Suj 2 = Nrj 2 . (9) 

There are N + 1 solutions of ft r , corresponding to the N + 1 normal collective modes. Let us for a moment suppress 
the index r of ft 2 . If ojq > 6ui 2 , all possible solutions for ft 2 are positive, physically meaning that the system oscillates 
harmonically in all its modes. If ui 2 , < Suj 2 , then a single negative solution exists. In order to prove this let us define 
the function 

N 1 

/(ft 2 ) = u 2 - scu 2 - ft 2 - v 2 n 2 ]T -3—02, (10) 

so that Eq. © becomes I (ft 2 ) = 0. We find that 



7(ft 2 )^oo as ft 2 -> -oo and 7(0) = lu 2 - 5lu 2 < 0, 

in the interval (— oo, 0]. As I (ft 2 ) is a monotonically decreasing function in this interval, we conclude that I (ft 2 ) = 
has a single negative solution in this case. This means that there is a mode whose amplitude grows or decays 
exponentially, so that no stationary configuration is allowed. Nevertheless, it should be remarked that in a different 
context, it is precisely this runaway solution that is related to the existence of a bound state in the Lee-Friedrichs 
model. This solution is considered in the framework of a model to describe qualitatively the existence of bound states 
in particle physics [18j |. 

Considering the situation where all normal modes are harmonic, which corresponds to the first case above, > Suj 2 , 
we define the renormalized frequency 

Oo 2 = uj 2 - Su) 2 = lim (w 2 - Nt] 2 ), (11) 

N— »oo 



in terms of which Eq. {5| in the limit N — > oo becomes, 

ft 

k=l 



In this limit, the above procedure is exactly the analog of the mass renormalization in quantum field theory: the 
addition of a counterterm — S^q 2 , allows one to compensate the infinity of to 2 , in such a way as to leave a finite, 
physically meaninful renormalized frequency Q. This simple renormalization scheme has been introduced earlier p^ |. 
Unless explicitly stated, the limit N — ► oo is understood in the following. 
Let us define a constant g, with dimension of frequency, by 



where Ato = ttc/R. The environment frequencies u>k are given by, 



(13) 
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"k = k— , fc=l,2,..., (14) 
where R is the radius of the cavity that contains the whole system. Then, using the identity 
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Eq. (|12p can be written in a closed form: 



:m\ ft c / ffi) 2 \ 

cot — = — + — 1 . (16) 

c / irg R\l \ irgc 



The solutions of the above equation with respect to ft give the spectrum of eigenfrequencies ft r corresponding to the 
collective normal modes. 
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In terms of the physically meaningful quantities £! r and ui, the transformation matrix elements turning the particle- 
field system to the principal axis are obtained. They are 



f — 



(r>2 - ^2)2 + £( 3r! 2 _ ^2) + ^2^2 



f 



(17) 



These matrix elements play a central role in the quantities describing the system. 

II. THE THERM ALIZATION PROCESS IN BARE COORDINATES 

We now consider the thermalization problem us ing bare coordinates. For the model described by Eq. (|TJ) this 
problem was addressed in an alternative way in [20j with the canonical Liouville-von Neumann formalism. We 
consider the initial state described by the density operator, 



p{t = 0) = po <g> pp 



(18) 



where po is the density operator of the particle, that in principle can be in a pure or in a mixed state and pp is the 
density operator of the thermal bath, at a temperature , that is, 



Pp = Z 1 exp 



-/3^w & (ala k + - 
fc=i ^ 



(19) 



with Zp = Hfcli z p being the partition function of the reservoir, and 



4 = Tr fe 



e -Puj k (ala k + l/2) 

The creation and annihilation operators given by 
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where uj^ — (<D,cjfe). The thermalization problem is addressed by investigating the time evolution of the state p{t). 

The thermalization problem concerns the time evolution of the initial state to thermal equilibrium. The subsystem 
corresponding to the particle oscillator is described by an arbitrary density operator po- As we will show, the 
expectation value of the number operator corresponding to particles will evolve in time to a value that is independent 
of the initial density operator poi the dependence will be exclusively on the mixed density operator corresponding to 
the thermal bath. 

Our aim is to obtain expressions for the time evolution of the expectation values for the occupation number and 
in particular for the one corresponding to particles. We will solve the problem in the framework of the Heisenberg 
picture. It is to be understood that when a quantity appears without the time argument it means that such quantity 
is evaluated at t = 0. The Heisenberg equation of motion for the annihilation operator a^(t) is given by 



d r - 

— a nit) = t \ H,a^(t) 

Due to the linear character of our problem, this equation is solved by writing a M (i) as 



(23) 



(24) 
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where all the time dependence is in the c -number functions B^ v (t). Then, Eq. (|23]) reduces to the following coupled 
equations for B^t): 

oo 

B^ (t) + Co 2 B^ (t) -Y / V^kB^ k (t) = 0, (25) 
fe=i 

oo 

B„ k (t) + ulB^t) - B^{t) Wk = 0. (26) 

fc=i 

These equations are formally identical to the classical equations of motion, Eqs. ((2]) and j3|), for the bare coordinates 
q^. Then we decouple Eqs. (f2"5|) and (f2"6")l with the same matrix {t r ^} that diagonalizes the Hamiltonian Eq. |T]). In 
an analogous manner, we write B^ v {t) as 

oo 

M^EW ( 27 ) 

such that from Eqs. (j2l)j) and (|2"6")) . we obtain the following equations for the normal-axis functions C r {t), 

c;(t) + nlc;(t) = o, (28) 

which gives the solution 

C;(t) = ^e^ + ^e- in ^. 
Then substituting this expression into Eq. ([?7|) we find 

oo 

= £ S ( a M ei ° rt + 6^-*°-*) . (29) 

The time independent coefficients a^, 6^ are determined by the initial conditions at t = Q for B^^t) and B^it). 
From Eqs. (f2"Tj) and (pM]) we find that these initial conditions are given by 



J tlV — /TT^> 



Buv — 



I CO 



~ V 2 

Using these equations, we obtain for and 

Co, 



B»u = \Hr 5 ^- ( 3 °) 



We write a M (t) and a* (t) in terms of a M and using Eqs. (f2T|) , ff22|) and (|24|) . 

OO 

= £ ( a ^W^ + P^(t)al) , (33) 

!/=0 

oo 

4w = EO^w^+^wsi) , (34) 
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where a^ v (t) and (3^, v (t) are the Bogoliubov coefficients given by, 



k-M ] -'n 

and 

Using the definition of B^^t) we get 

oo ; j.r j.r 



a^{t) = - 1 =B llu {t) - iJ-^B^t) (35) 



/V W = - 1 =B^(t) + iJ-^B^t) . (36) 



/ f' r f r ( o 



r=0 

+ [(Q r - ^)e in ^ + (Q r + wpje-^*] } , (37) 



and 



00 rri—+r+r i w 



- [(Sir - w^e* 1 -* + (Q r + u^e^*] } . (38) 
Now we study the time evolution of n^(t), the expectation value of the number operator N^t) — at(t)a^(t), that is, 

n^t) = Tr[al(t)a^t)p ®p p ] . (39) 
Using the basis |no, ni, n2, ...njv) we obtain, 



X 

u=0 u=0 



.(*) = E [i«^wi a + i^wi 2 ] ^ + E i/w*)i a - (4°) 

where 



oo 



n = / , n MPoN ( 41 ) 

n=0 

is the expectation value of the number operator corresponding to the particle and the set {n k } stands for the thermal 
expectation values corresponding to the thermal bath oscillators, given by the Bose-Einstein distribution, 

n k = -z— • (42) 

In Eq. ([40]) there appears a term that does not depend on the temperature of the thermal bath. This term has its 
origin in the instability of the initial bare vacuum state, |0, 0, 0). To see this, we compute the expectation value of 
the time dependent number operator N^t) = at (t)a^(f) in this vacuum state. Thus all the terms containing operators 
different from the identity give a zero contribution. The only term, that gives a non-zero contribution comes from the 
normal ordering and is just the last one in Eq. (|40|) . This term leads to the creation of excited states (particles, in a 
field theoretical language) from the initial unstable bare vacuum state. 

We are interested in evaluating the expectation value of the number operator corresponding to particles. Thus 
taking fi — in Eq. (|4"0)) and using Eq. (|4*2"|) , we obtain 



no(t) = [\aoo(t)\ 2 + \Mt)\ 2 ]n + J2[\a ok (t)\ 2 + \/3 Qk (t)\' 



k=l 



J3u k _ i 



-\0oo(t)\ 2 +J2\P 0k (t)\ 2 , (43) 



fc=i 
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where the coefficients of this expression are [20j , 

e -*gt/2 



and 



such that 



"oo (t) 

#»(*) 

otok(t) 



16lok 



(2u) + 2k — ing) e 



2 —%Kt 



2 Ant 



(2u> — 2k — iirg) e 



nge 



-irgt/2 



[{■Kg + 2k) e- lKt - {irg - 2k) e lKt ] , 



8ujK 

[jj^ (g + uj k )^/g^J e~ tuJkt [Z^s/2g~Eu 
2uj (oj? — uj 2 + iirgojk) 
{2k + 2uj- ing) ^ iKt 
{2k — 2iUk — mg) 



uj 4k 
{2u> -2k- ing) ^ Kt 
{2k + 2u>k + iirg) 



-ngt/2 



ui k {ujk - uj)^fgKde lbJkt fu* ' ^/2g~Kd 



2u {to 2 — up- — ingujk) 
(2w + 2k - ing) lKt 



uj 4k 



(2k + 2ujk — i^g) 



{2uj — 2k — ing) 

- -i 

(2k — 2ujk + ing) 



-ngt/2 



(44) 
(45) 

(46) 



(47) 



k = a/w 2 - 7r 2 g 2 /4. 



(48) 



The parameter k measures the intensity of the interaction: if k 2 >> 0, i.e g << 2uj /n, we are in the weak coupling 
regime. On the contrary if k 2 « 0, i.e g » 2uj/t:, the system is in the strong coupling regime. Here we will restrict 
ourselves to the weak coupling regime. This case includes the important class of electromagnetic interactions, g — auj, 
with a being the fine structure constant a = 1/137 [T3 |. 



In the continuum limit Auj 

n {t), 



0, sums over k become integrations over a continuous variable uj and we obtain for 



where 



n {t) 



-ngt 



-4 + ^_9_ ^ 2 -2 _ ^2^ cos (2 K £) - - 9 K sin(2^) 



no 



Tr 2 g 2 e-^ t 
V6uj 2 k 2 



[2d 2 + {2uj 2 - TT 2 g 2 ) cos{2Kt) - 27rgKsin(2frf)] 



F{uj,uj,g,t) 
(e^ - 1) 



G{uJ,u,g,t) 



F(uj,uj,g,t) 



uj(u? +Q 2 ) 



[{uj 2 — UJ 2 ) 2 + TT 2 g 2 UJ 2 

{uj 2 -uj 2 ) 



-wgt 



4k 2 



\4uj 2 - n 2 g 2 cos(2fi-i) 



-2ngK 



{uj 2 + UJ 2 ) 
4ujuj 2 



sin(2«t)] — 



-irgt/2 



- [2k cos(wi) cos{Kt) 



{UJ 2 + UJ 2 ) 



sm{ujt) sin(/d;) —ng 



{uj 2 - uj 2 ) 

{UJ 2 +LU 2 ) 



cos{ujt) sm{Kt)] 



(49) 



(50) 



and 



G{uj,uj,g,t) = 



uj{uj — u>Y 



[{uj 2 - UJ 

{uj 2 +uj 2 ) 



1 



,-irgt 



4k 2 



4uj 1 



2n 2 g 2 ujuj 
{uj - uj) 2 



2 2 1 
-7T g 



[UJ - UJ 



cos{2Kt) — 2ircjK- — _ sin(2/«i) 
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e -fgt/2 

[2k cos(uit) cos(nt) — 2u> sintuit) sin(/rf) 

K 

(lU + Cu) 1 

— 7T(7- — cos(wi) sm(nt)\ > . (51) 

It is to be noted that the second and the third lines in Eq. ([49]) are independent of the initial distributions. Also the 
integral in the third line of G(w, u>, g, t) is logarithmically divergent. We can understand the origin of this divergence in 
the following way: suppose that initially, in the absence of the linear interaction, we prepare the system in its ground 
state, that is, at t = we have |0,0, ...,0). Then, we can compute, in the Heisenberg picture, the time evolution 
for the expectation value of the number operator corresponding to the particle, that is (0, 0, 0|ctj(£)ao(£)|0, 0, 0). 
Taking /i = we obtain, 



(0, 0, ...,0|a$(i)fio(i)|0, 0, 0) = |/3 00 (i)| 2 + £ |A>k(*)l a , ( 52 ) 

k=l 

which in the continuum limit gives the second line of Eq. (|49[) . Then, the origin of the divergence appearing in Eq. (|49p 
is interpreted as the excitations produced from the unstable bare (vacuum) ground state, as a response to the linear 
interaction. 

As we are interested in the thermal behavior of no(t) only, the second line and the term G(u>,u>, g,t) in the third 
line of Eq. (|49|1 can be neglected. This is a renormalization procedure. Thus we write the following renormalized 
expectation value for the particle number operator, 



n (t) = K(u,g,t)n - 



F(lu, uj, g, t) 
(eP» - 1) 



(53) 



where 



K(u,g,t) 



^-f - (2Q 2 - ttV) cos(2frf) - ^ 9 K sin(2«i) 



(54) 



In the limit t — > oo, fio(t) has a well defined value, that is, the system reaches a final equilibrium state. Also, 
since K(u>, g,t — > oo) — > 0, this final equilibrium state is independent of no. The equilibrium expectation value of the 
number operator corresponding to the particle is independent of its initial value, and the only dependence is on the 
initial distribution of the thermal bath, that is, the particle thermalizes with the environment. Before the interaction 
enters into play for t < 0, nit < 0) = tiq, then we have that K(ui,u),g,t < 0) = 1. Taking t = in Eq. (fB"4"|) we 
obtain K (u>, u>,g, t = 0) = Co 2 / k 2 + ir 2 g 2 {2uj 2 — it 2 g 2 ) / (8uj 2 n 2 ) . Thus K(ui, u, g, t) is a discontinuous function of t; the 
discontinuity appearing just at t = 0. From the physical standpoint this discontinuity can be viewed as a response to 
the sudden onset of the interaction between particles and the environment. 

Although the integral in Eq. (|53[) can not be computed analytically, we can perform numerical calculations, for 
example in Fig. [T]we display the time behavior for no = 1, u> = 1, j3 = 2 and g — 0.1; (i > 1). In the next section 
we develop an alternative approach based on the notion of dressed particles. We will find that, in this new realm, no 
renormalization is needed. 



III. DRESSED COORDINATES AND DRESSED STATES 



Let us start with the eigenstates of our system, |no, n\, ri2---) , represented by the normalized eigcnfunctions in terms 
of the normal coordinates {Q r }, 



un a ...(Q,t) = n 




H„ 



r e" 



where H Ug stands for the n s -th Hermite polynomial and r is the normalized vacuum eigenfunction, 

1 o = Ae 2 ^r=o r ^ r 



(55) 



(56) 



5 10 15 20 25 30 35 40 

I 

FIG. 1: Time behavior for n (t) given by Eq. (J53J) for (t > 1), no = 1, O = 1, f3 = 2 and g = 0.1 



We introduce dressed or renormalized coordinates q' and {q^} for, respectively, the dressed particle and the dressed 
field, defined by, 

=XX"^V2r, (57) 
r 

valid for arbitrary R and where <D M = {(D, w^}. In terms of dressed coordinates, we define for a fixed instant, t = 0, 
dressed states, \ko,k\, k-2...) by means of the complete orthonormal set of functions 

ipKoKx.-Xq') = II 

where q'^ — {q' , <^}, — {iD, cjj}. Notice that the ground state T in the above equation is the same as in Eq. (|55[) . 
The invariance of the ground state is due to our definition of dressed coordinates given by Eq. (fS"T|) . Each function 
V'reorei. ••(#') describes a state in which the dressed oscillator q'^ is in its K M -th excited state. 

It is worthwhile to note that our renormalized coordinates are new objects, different from both the bare coordinates, 
q, and the normal coordinates Q. In particular, the renormalized coordinates and dressed states, although both are 
collective objects, should not be confused with the normal coordinates Q, and the eigenstates Eq. ([53)1 . While the 
eigenstates <fi are stable, the dressed states ip are all unstable, except for the ground state obtained by setting {k m = 0} 
in Eq. ([55)) . The idea is that the dressed states are physically meaningful states. This can be seen as an analog of 
the wave-function renormalization in quantum field theory, which justifies the denomination of renormalized to the 
new coordinates q' . Thus, the dressed state given by Eq. (|58|) describes the particle in its Ko-th excited level and 
each mode k of the cavity in the — th excited level. It should be noticed that the introduction of the renormalized 
coordinates guarantees the stability of the dressed vacuum state, since by definition it is identical to the ground state 
of the system. The fact that the definition given by Eq. (|57p assures this requirement can be easily seen by replacing 
Eq. (|57|) in Eq. (|58|) . We obtain Tv{q') oc ro(Q), which shows that the dressed vacuum state given by Eq. (|58|) is the 
same ground state of the interacting Hamiltonian given by Eq. ((5|) . 

The necessity of introducing renormalized coordinates can be understood by considering what would happen if we 
write Eq. (|58|) in terms of the bare coordinates q^. In the absence of interaction, the bare states are stable since 
they are eigenfuntions of the free Hamiltonian. But when we consider the interaction they all become unstable. The 
excited states are unstable, since we know this from experiment. On the other hand, we also know from experiment 
that the particle in its ground state is stable, in contradiction with what our simplified model for the system describes 
in terms of the bare coordinates. So, if we wish to have a nonperturbative approach in terms of our simplified 
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model something should be modified in order to remedy this problem. The solution is just the introduction of the 
renormalized coordinates q' as the physically meaningful ones. 

In terms of bare coordinates, the dressed coordinates are expressed as 



where 



If we consider an arbitrarily large cavity (i? — ► oo), the dressed coordinates reduce to 



ifx, = X) "V?"' ( 59 ) 



=i;^v^. (so) 



with A 00 (u,g) given by, 



la = A oo{u, g)q , (61) 

q[ = a, (62) 



A 00 (^g)-^=J^ (fi2 _- 2)2+ ^ w - (63) 



In other words, in the limit i? — > oo, the particle is still dressed by the field, while for the field there remain bare 
modes. 

Let us consider a particular dressed state |r^(0)), represented by the wavefunction ?/W--i(ai)o---(<z')- ^ describes the 
configuration in which only the dressed oscillator q'^ is in the first excited level. Then the following expression for its 
time evolution is valid fl3| : 

Kit)) = E^(*)i r i(°)) 

v 

s 

Moreover we find that 

J2\F V (t)f = l- (65) 

V 

Then the coefficients /'""(t) are simply interpreted as probability amplitudes. 

In approaching the thermalization process in this framework, we have to write the initial physical state in terms of 
dressed coordinates, or equivalently in terms of dressed annihilation and creation operators a 1 ^ and ajj instead of 
and ajj. This means that the initial dressed density operator corresponding to the thermal bath is given by 



Pa = z 8 1 ex P 



^E Wfc \ a 'k a 'k + 



k=l 



(66) 



where we define 




Now we analyze the time evolution of dressed coordinates. 
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IV. THERMAL BEHAVIOR FOR A CAVITY OF ARBITRARY SIZE WITH DRESSED COORDINATES 



The solution for the time-dependent annihilation and creation dressed operators follows similar steps as for the bare 
operators. The time evolution of the annihilation operator is given by, 



d 
dt 



and a similar equation for a'Ht). We solve this equation with the initial condition at t = 0, 



(69) 



<(0) 



(70) 



which, in terms of bare coordinates, becomes 



N 



a U°)= E vf^+ 



r,u— 



"t- ...... - lb [i L i/ 



(71) 



We assume a solution for a' (t) of the type 



u=0 



(72) 



Using Eq.flJ we find, 



iCl r t\ 



(73) 



r=0 



In the present case the time independent coefficients are different from those in the bare coordinate approach, Eq. (|29p . 
The initial conditions for B' (t) and B' (t) are obtained by setting t — in Eq. f72j) and comparing with Eq. (fTTj) : 



Then 



5^(0) = 



r=0 



r=0 



(74) 
(75) 



Using these initial conditions and the orthonormality of the matrix {t^} we obtain = 0, b'£ = it r fl /y / 2fl r . Replacing 
these values for a'l and b'£ in Eq. ([73]) we get 



OO 17. ,7. 

Vv -iO r t 



r=0 



(76) 



We have 



iV 



»;.(*) = E 

r,f — 
iV 



/2SI 



r,v= 



-in r t 



E-^(*K > 



(77) 



where 



i(l r t 



(78) 
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For the occupation number n'^(t) = (a'J (t)a' (£)) we get 

n' fl (t)=Tr(a'}(t)a' fl (t)p' o ®p' ) . (79) 

where p' is the density operator for the dressed particle and p'^ is the density operator for the thermal bath, which 
coincides with the corresponding operator for the bare thermal bath if the system is in free space (in the sense of an 
arbitrarily large cavity) [H, • 

To evaluate n'^t) we choose the basis |no,ni, ...,tiat) = Jl^Lo l n /j)' where \n^) are the eigenvectors of the number 
operators ajja^. From Eq. ([77)) we get 

oo 

v,p=Q 
oo 

= \Uv(t)\ 2 &ya'„ + e f; p {t)u v {t)a'y a ' p . m 

In the basis |no,rti,rt2, • • •) we obtain, 



^(*) = i/^(t)ix+Ei^(*)i v *' ( 81 ) 

k=l 

where n' and n' k are the expectation values of the initial number operators, respectively, for the dressed particle and 
dressed bath modes. We assume that, dressed field modes obey a Bose-Einstein distribution. This can be justified 
by remembering that in the free space limit, R — > oo, dressed field modes are identical to the bare ones, according 
to Eqs. (|61[) and (|62|) . Now, no term independent of the temperature appears in the thermal bath. This should be 
expected since the dressed vacuum is stable, particle production from the vacuum is not possible. Setting p = in 
Eq. (|8ip we obtain the time evolution for the ocupation number of the particle, 

oo 

^(*)H/oo(*)lX + £l/ofc(*)|V fe . (82) 
fe=i 

V. THE LIMIT OF ARBITRARILY LARGE CAVITY: UNBOUNDED SPACE 

In a large cavity (free space) we must compute the quantities /oo (t) and fok (t) in the continuum limit to study the 
time evolution of the ocupation number for the particle. Remember that in Eqs. (|17[) . LOk = knc/R, k = 1,2, ... and 
r] = \/2gSu), with Alu = (u>i+i — LOi) = ixcjR. When R — > oo, we have Aw — > and AQ — > and then, the sum in 
Eq. (|78p becomes an integral. To calculate the quantities f^, v {t) we first note that, in the continuum limit, Eq. (|17p 
becomes 

t5 - $VMl = lim Vly/2gATL 

S - *£§lgVSB. (84) 

In the following, we suppress the labels in the frequencies, since they are continuous quantities. 
We start by defining a function W(z), 

oo 2 2 

IF(z) = z 2 -^ + £4^. (85) 

We find that the fi's are the roots of W(z). Using rf = 2gAa;, we have in the continuum limit, 

W(z) = z 2 -0j 2 + 2gz 2 f (86) 
Jo - z l 
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For complex values of z the above integral is well denned and is evaluated by using Cauchy theorem, to be 



W{z) 



_ J z 2 + ignz — lu 2 , Im(z) > 
z 2 — igirz — lu 2 , Im < . 



We now compute /oo(i) = X)r=o(^o) e which, in the continuum limit, is given by 



/oo(t) 



(#) 2 e~ lUt dn 



We find that, 



1 



W(Q) ' 

and since the fi's are the roots of W(z), we write Eq. I|88p as 

, , , If dze~ lzt 

17T/ C VF(z) 



(87) 



(88) 



(89) 



(90) 



where C is a counterclockwise contour in the z-plane that encircles the real positive roots of W(z). Choosing a contour 
infinitcsimally close to the positive real axis, that is z = a — ie below it and z = a + ie above it with a > and 
e — ► + , we obtain 



1 [°° 
/oo(*) = — daae~ 
™ Jo 



1 



1 



W(a-ie) VF(a + ie) 



In the limit e — > + , Eq. ([57]) gives VF(a ± ie) = a 2 — cl> 2 ± igira which leads to 

/oo(*) = Ci{t;ui,g) + iS^Q, g), 

where 

a 2 cos(crf) 
(a 2 -co 2 ) 2 + it 2 g 2 a 2 ' 



f OO 

Ci{t;w,g) = 2g da- 
Jo (" 

Si(t;uj,g) = -2g 



00 , a 2 sin(ai) 
aa- 



(a 2 — lu 2 ) 2 + ir 2 g 2 a 2 



(91) 
(92) 

(93) 
(94) 



Notice that C\(t = 0;<D,g) = 1 and S\(t = 0;0u,g) = 0, so that foo(t = 0) = 1 as expected from the orthonormality 
of the matrix (i^). The real part of foo(t) is calculated using the residue theorem. For k 2 = Co 2 — ir 2 g 2 /A > 0, which 
includes the weak coupling regime, one finds 



Cx(t;uj,g) = e-* gt l 2 [cos(/si) - — sux(/ct)l (k 2 > 0). 

L 2k J 

Although Si(t;cj,g) cannot be analytically evaluated for all t, however for long times, i.e. t 1/u), we have 



Si(t;u,g) 



4 9 , is 



Thus, we get for large t 



1/ooWI 



2 ^ -irnt 



7TO . , . 
cos{K,t) — — sm(Kt) 



+ 



16g 2 



Next we compute the quantity /ofc(i) — Yl r =o m ^ ne continuum limit. It is 



{lu 2 - ft 2 ) ~ m Jc {lu 2 - z 2 )W{z) ' 
where r\ = \f2gSuj. Taking the same contour as that used to calculate foo{t), we obtain 



ITT 



da 



W{a - ie)[{a - ie) 2 - lu 2 ] W{a + ie)[(a + ie) 2 - lu 2 ] 



(95) 



(96) 



(97) 



(98) 



(99) 
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FIG. 2: Time behavior for n' (t) given by Eq. for (i > 1), no = 1, Q = 1, f) = 2 and g = 0.1 



Thus, taking e — > + fou(t) is written as 

/ou(f) = wVAcJ [C 2 (uj,t;uj,g) + iS 2 (uj,t;Cu, g)] 

where 

a 2 cos(at) 



C 2 (u,t;Q,g) = (2g)v 
S 2 (uJ,t;u},g) = -(2g) 



da- 



o ' (uj 2 - a 2 ) [(a 2 - ui 2 ) 2 + ir 2 g 2 a 2 ] 
a 2 sin(ai) 



da 



q (w 2 — a 2 ) [(a 2 — cj 2 ) 2 + 7r 2 <7 2 a 2 ] 



(100) 

(101) 
(102) 



Notice that the integrals defining the functions C 2 and S 2 are actually Cauchy principal values. 
The function C 2 is calculated analytically using Cauchy theorem; we find 



C 2 (oJ,t;u,g) 



2g 
ng 



,-Tvgt/2 



LO 2 + uj 2 



2k (uj 2 - uj 2 ) 2 + Ti 2 g 2 uj 2 
nguj 



{uj 2 — uj 2 ) 2 + n 2 g 2 uj 2 
sin Kt 



cos Kt 



sin ujt 



(103) 



(uj 2 — UJ 2 ) 2 + TT 2 g 2 LJ 2 

The function S 2 cannot be evaluated analytically for all t, it has to be calculated numerically. For long times, we have 



S 2 (t;uj,g)^ / (*>-)• 

UJ z UJ^t a UJ 



(104) 



In the continuum limit, we get the average of the particle ocupation number, 

poo 

n' (t) = [Cj(t;uj 7 g) + S 2 1 (t;uj 7 g)]n' Q + / <W [C 2 2 (w,t; Q,g) + %(u,t;Q,g)] n'(cv) , (105) 



where n'(uj) = l/(e* w — 1) is the density of occupation of the environment modes, the functions C\ and C 2 are given 
by Eqs. (|95|) and (|103[) while the functions S\ and are given by the integrals Eqs. (|94|) and (|102[) . respectively. In 
Fig. [2] we display the behavior in time for n = 1, Q = 1, /? = 2 and g = 0.1; (t > 1). 

The important point, that is seen from Fig.Q]and Figgis that, for long times, both the bare and dressed ocupation 
numbers of the particle approach smoothly to the same asymptotic value. Moreover this value is the one expected on 
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physical grounds, obtained from the Bose distribution at the final equilibrium temperature. In fact, taking [3 — 2 and 
u> = 1, as used in the plots, we get 

riociu)) = l/(e^ - 1) = 0.156 

Therefore both methods, and in particular our dressed state formalism describes very precisely the thermalization 
process. 

VI. FINAL REMARKS 

We have considered a linearized version of a particle-environment system and we have carried out a non-perturbative 
treatment of the thermalization process. We have adopted the point of view of renouncing to an approach very close to 
the real behavior of a nonlinear system, to study instead a linear model. As a counterpart, an exact solution has been 
possible. This realises a good compromise between physical reality and mathematical reliability. We have presented 
an ohmic quantum system consisting of a particle, in the larger sense of a material body, an atom or a Brownian 
particle coupled to an environment modelled by non-interacting oscillators. We have used the formalism of dressed 
states to perform a non-perturbative study of the time evolution of the system, contained in a cavity or in free space. 
Distinctly to what happens in the bare coordinate approach, in the dressed coordinate approach no renormalization 
procedure is needed. Our renormalized coordinates contain in themselves the renormalization aspects. As far as the 
thermalization process is concerned from a physical viewpoint, both bare and dressed approaches are in agreement 
with what we expect for this process. For long times, all the information about the particle occupation numbers 
depends only on the environment. Both curves in Fig. [T] and Fig. [2] approach steadly to an asymptotic value of the 
bare and dressed ocupation numbers of the particle, which is the physically expected one at the given temperature. 
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